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THE SHARP INTERFACE LIMIT FOR THE 
STOCHASTIC CAHN-HILLIARD EQUATION 

D.C. ANTONOPOULOU®*, D. BLOMKER*, AND G.D. KARALI f * 


Abstract. We study the e-dependent two and three dimensional stochastic Cahn-Hilliard equation 
in the sharp interface limit e —> 0. The parameter e is positive and measures the width of transition 
layers generated during phase separation. We also couple the noise strength to this parameter. 
Using formal asymptotic expansions, we identify the limit. In the right scaling, our results indicate 
that the stochastic Cahn-Hilliard equation converge to a Hele-Shaw problem with stochastic forcing 
on the curvature equation. In the case when the noise is sufficiently small, we rigorously prove that 
the limit is a deterministic Hele-Shaw problem. Finally, we discuss which estimates are necessary 
in order to extend the rigorous result to larger noise strength. 

Resume: Nous etudions 1’ equation Cahn-Hilliard stochastique dependante en e, posee en dimen¬ 
sions deux et trois dans la limite de l’interface nette e —> 0. Le parametre e est positif et mesure la 
largeur de couches de transition generees pendant la separation de phase. Nous avons aussi couple 
la puissance de bruit a ce parametre. A l’aide de series asymptotiques formelles, nous determinons 
la limite. Dans l’echelle propre, nos resultats indiquent que 1’ equation Cahn-Hilliard stochas¬ 
tique converge vers un probleme Hele-Shaw avec une force stochastique presente a l’equation de la 
courbure. En cas de bruit suffisamment petit, nous prouvons rigoureusement que la limite est un 
probleme Hele-Shaw deterministe. Finalement, nous discutons quelles estimations sont necessaires 
ahn de prolonger le resultat rigoureux en presence de bruit d’une puissance plus grande. 


1. Introduction 


In this paper we consider the sharp interface limit of the stochastic Cahn-Hilliard equation 
(1.1) d t u = A(— eAu + e _1 / / (tt)) + e a W(x,t ), 


for times t 

d G {2,3} 

( 1 . 2 ) 


G [0,T] subject to 


Neumann boundary conditions on a bounded domain in V C M. d , 


du 

dn 


dAu 

dn 


on dV . 


Here u : T> x [0, T] —> R is the scalar concentration field of one of the components in a separation 
process, for example of binary alloys. 
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We always assume that the spatial domain T> has a sufficiently piece-wise smooth boundary. The 
typical example for the bistable nonlinearity is f'(u) := d u f{u) = u(u 2 — 1), where the primitive 
/ := \{u 2 — l) 2 is a double well potential with equally deep wells taking its global minimum value 
0 at the points u = ±1. The small parameter e measures an atomistic interaction length that fixes 
the length-scale of transition layers between 1 and —1. Obviously, the solution u depends on e, 
which we usually suppress in the notation, but wherever needed we denote the solution u by u £ . 

The forcing W denotes a Gaussian noise, which is white in time and possibly colored in space. 
Finally, the noise-strength e a scales with the atomistic length scale, where o is some parameter. 

In applications, there are often different sources of the noise. If we split the noise e a W = A V+W, 
then the term V arises from fluctuations in the chemical potential, while W models fluctuations 
directly in the concentration. Here, V needs to be a spatially smooth space-time noise since it is 
under the Laplacian, while the additive noise W might be rougher. 

The stochastic Cahn-Hilliard equation is a model for the non-equilibrium dynamics of metastable 
states in phase transitions, [261 [32 , 35] . The deterministic Cahn-Hilliard equation with W = 0 
has been extended to a stochastic version by Cook, {26J (see also in [35]), incorporating thermal 
fluctuations in the form of an additive noise. Such a generalized Cahn-Hilliard model, [36], is based 
on the balance law for micro-forces. In this case, the additive term V in the chemical potential is 
given by fluctuations in an external field. See |421136] : cf. also [33], where the external gravity field 
is modeled. The noise W, which is independent of the free energy, stands for the Gaussian noise 
in Model B of [32], in accordance with the original Cahn-Hilliard-Cook model. 

We can model the noise either as the formal derivative of a Wiener process in the sense of Da 
Prato &; Zabczyck [29] or as a derivative of a Brownian sheet in the sense of Walsh, [50]. Dalang & 
Quer-Sardanyons m showed that both approaches are actually equivalent. Thus, we focus on the 
Q-Wiener-process in the sense of [29] . which may be defined by a series in the right orthonormal 
basis with independent coefficients given by a sequence of real valued Brownian motions. 

Existence of stochastic solution for the problem (11.11) has been established under various as¬ 
sumptions on the noise-terms, see for example in [27l 28]: we also refer to the results presented in 

[ZI] EZl EE HB]. 

1.1. Phase transitions and noise. Concerning the Cahn-Hilliard equation posed in one dimen¬ 
sion, in [8] the authors analyzed the stochastic dynamics of the front motion of one dimensional 
interfaces. In the absence of noise we refer also to 01 for the dynamics of interfaces and the 
construction of a finite dimensional manifold parametrized by the interface positions. This is a key 
tool for studying the stochastic case, but it fails in dimension two and three, as interfaces are no 
longer points, but curves or surfaces respectively. 

An interesting result is m, where, on the unbounded domain, a single interface moves according 
to a fractional Brownian motion, which is in contrast to the usual Brownian motion in most of the 
other examples. Note that the one-dimensional case is significantly simpler, since the solutions 
can be fully parametrized by their finitely many zeros. Therefore, one needs only to consider the 
motion of solutions along a finite dimensional slow manifold, and the stability properties along 
such a manifold; see for example m- Here we try to follow similar ideas, despite the fact that 
the driving manifold is infinite dimensional, and parametrized by the curves or surfaces of the zero 
level set. 

A slightly simpler model, due to the absence of mass-conservation, is the stochastic Allen Cahn- 
equation, the so called Model A of m • The interface motion of stochastic systems of Allen-Cahn 
type have been analyzed in m- In [35[ 19] . the authors studied the stochastic one-dimensional 
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Allen-Cahn equation with initial data close to one instanton or interface and proved that, under 
an appropriate scaling, the solution will stay close to the instanton shape, while the random per¬ 
turbation will create a dynamic motion for this single interface. This is observed on a much faster 
time scale than in the deterministic case. This result has been also studied in [ 52] via an invariant 
measure approach. 

If the initial data involves more than one interfaces, it is believed that these interfaces exhibit 
also a random movement, which is much quicker than in the deterministic case, while different 
interfaces should annihilate when they meet, [32]. We also refer to [49] or [37]. The limiting 
process should be related to a Brownian one (cf. in [33] for formal arguments). A full description 
of all the ideas for the analysis of the interface motion based on [25 ] and [|8j is presented in m- 
In [48] the authors considered the stochastic Allen-Cahn equation driven by a multiplicative noise; 
they prove tightness of solutions for the sharp interface limit problem, and show convergence to 
phase-indicator functions; cf. also in [52] for the one-dimensional case with an additive space-time 
white noise for the proof of an exponential convergence towards a curve of minimizers of the energy. 

The space-time white noise driven Allen-Cahn equation is known to be ill-posed in space dimen¬ 
sions greater than one, [50i 129] . and a renormalization is necessary to properly define the solutions. 
We refer to [38] 00] for more details. 

For a multi-dimensional stochastic Allen-Cahn equation driven by a mollified noise, in [39] . it 
is shown that as the mollifier is removed, the solutions converge weakly to zero, independently of 
the initial condition. If the the noise strength converges to zero at a sufficiently fast rate, then 
the solutions converge to those of the deterministic equation, and the behavior is well described 
by the Freidlin-Wentzell theory. See also [T6]. In m, for a regularized noise, which is white in 
the limit, the author extending the classical result of Funaki |33j to spatial dimensions more than 
two, derived motion by mean curvature with an additional stochastic forcing for the sharp interface 
limit problem. Recently, in [T], for the case of an additive ‘mild’ noise in the sense of [Ml El], 
the first rigorous result on the generation of interface for the stochastic Allen-Cahn equation has 
been derived; the authors proved layer formation for general initial data, and established that the 
solution’s profile near the interface remains close to that of a (squeezed) traveling wave, which 
means that a spatially uniform noise does not destroy this profile. 

1.2. Main Results. Due to phase separation, the solution u of the stochastic Cahn-Hilliard equa¬ 
tion m, which is related to the mass concentration, tends to split in regions where u ~ ±1 
and with inner layers of order 0(e) between them. We shall study the motion of such layers in 
their sharp interface limit e —> 0. The rigorous complete description of the motion of interfaces 
in dimensions two and three stands for many years as a wide open question. With this paper, we 
tried to contribute towards a full answer by providing first some formal asymptotics for the case of 
general noise strength. We were able to identify the intermediate regime a = 1 when in the limit 
e —> 0, the limiting solutions satisfy a stochastic Hele-Shaw problem, (see (12.51) ]). While for smaller 
noise strength a > 1 we conjecture that the limit solves a deterministic Hele-Shaw problem. 

In addition, as a first rigorous step, we derive the sharp interface limit given by the deterministic 
Hele-Shaw problem in the case when the noise is sufficiently small, (see Theorem 13.101) . Unfortu¬ 
nately, we are not able to treat the case <r close to 1, we refer to Remark 13.111 A key problem 
is that spectral problems of the linearized operators are not yet understood in the general setting 
which is necessary for the proof, (see Remark 13.131) . 

The result for the Cahn-Hilliard equation in the deterministic case (when W = 0) has been 
already studied in [2], Given a solution of the deterministic Hele-Shaw problem, the authors 
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constructed an approximation of (ll.ljl without noise admitting this solution, as the interface moves 
between ±1. The analysis thereof will be the foundation of some of our results. The main technical 
problem of this strategy, is that the manifold of possible approximations is parametrized by an 
infinite dimensional space of closed curves or surfaces. Furthermore, the spectral information 
provided so far for the linearized problem, necessary for a qualitative study of the approximation is 
insufficient. This is due mainly to the fact that most of the larger eigenvalues actually are related 
to the fast motion of the interface itself, which do not obstruct the stability of the manifold. Thus, 
this approximation can only be valid on time-scales of order 0(1). We refer to the main result in 
Theorem ITlOl where we approximate solutions of the stochastic Cahn-Hilliard equation by solutions 
of the deterministic Hele-Shaw problem only for very small noise and only on ©(^-time-scales. 

A simpler case is the motion of droplets for the two-dimensional Cahn-Hilliard or the mass 
conserving Allen-Cahn equation. Here, the solutions can be fully parametrized by finite dimensional 
data, namely the position and radius of the droplets. See mm, and [7] (stochastic problem), for 
droplets on the boundary, and mm for droplets in the domain in the absence of noise. The 
approximation in these cases is valid for very long time-scales. 

1.3. Outline of the paper. In Section [2l we present the formal derivation of the stochastic Hele- 
Shaw problem from and identify the order of the noise strength that leads to a nontrivial 
stochastic modification of the deterministic limiting problem. 

In Section [3l we provide a rigorous definition of the setting and state the main results, which 
we then prove in Section [fj We concentrate on small noise strength and show in that case, that 
solutions of (tHD in the sharp interface limit of £ —>• 0 are well approximated by a Hele-Shaw 
problem. We will see that the main limitation towards a better approximation result is the lack of 
good bounds on the linearized operator. 


2. Formal asymptotics 


In this section we present some formal matching asymptotics applied to m that establish a 
first intuition towards a rigorous proof for the stochastic sharp interface limit. We remark first that 
the stochastic C-H equation (11.11) can be written as a system, where v is the chemical potential. 


( 2 . 1 ) 


dt.u = —Av + £ a yV, 



/'(«) 


+ eAu. 


Formally in m, and later more rigorously in E3 using a Hilbert expansions method, the as¬ 
ymptotic behavior for e -A 0 of the following deterministic system has been analyzed: 


d t u £ = —Av £ + G\, 

v £ = -1 ^aL +£ Au £ + G 2 , 

£ 

where now G\(x,t\e) and G 2 (x,t\£) are deterministic forcing terms. The sharp interface limit 
problem in the multidimensional case demonstrated a local influence in phase transitions of forcing 
terms that stem from the chemical potential, while free energy independent terms act on the rest 
of the domain. In addition, the forcing may slow down the equilibrium. 
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Given an initial smooth closed n— 1 dimensional hypersurface To in T> (this definition covers 
also the union of closed interfaces) then the limiting chemical potential 

(2.2) v := lim (eA u £ - + G 2 ), 

£—>■0+ 

satisfies the following Hele-Shaw free boundary problem (assuming that the limits exist) 

' Av = lim G\ in V\T(t), t > 0, 

£->0+ 

d n v = 0 on dV, 

on v = A H + lim G 2 on T(f), 

£—>0+ 

V = 7 j( d n v + - d n v~) on F(t), 

k r(o) = r 0 , 

where T(t) is the zero level surface of the limiting u(t), which is for fixed time t a closed n— 1 
dimensional hypersurface of mean curvature H = H(t) and of velocity V = V(t) that divides the 
domain T> in two open sets T> + (t ) and D~(t). The constant A is positive, and n is the unit outward 
normal vector at the inner and outer boundaries. 

According to the aforementioned arguments, each perturbation Gi has a different physical mean¬ 
ing and appears in a different equation when C-H is presented as a system. We will use some of 
the ideas of the deterministic asymptotic analysis, but we will see in the following that when the 
terms Gi are noise terms and small in e, they still have an impact on the limiting behavior. 


2.1. Formal derivation of the stochastic Hele-Shaw problem. In order to observe the limit 
behavior of (11.ip at larger noise strengths, we fix a = 1 as we expect the noise strength e to be the 
critical one. In order to avoid calculations with formal noise terms, where the order of magnitude 
and the definition of products is not always obvious, we use the following change of variables 

u £ := u £ - eW , 


where we assume that the Wiener process is spatially sufficiently smooth. Taking differentials, it 
follows that u £ and v £ solve the system 


(2.4) 


d t u £ = -Av e , 

\ 

v £ = — f'(u £ + eW) + eAu £ + £ 2 AW. 
£ 


Observe that for spatially smooth noise on time-scales of order 1 


u £ = u £ + 0(e), 

is an approximate solution of m for small e. Furthermore, the main advantage of the above 
representation as a system is based on the fact that (12.41) is now a random PDE without stochastic 
differentials and all terms appearing are spatially smooth and in time at least Holder-continuous. 
Thus, we can treat all appearing quantities as functions and analyze the equation path-wisely, i.e., 
for every fixed realization of the underlying Wiener process W. 
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Therefore, we are able to follow the ideas of the formal derivation presented in m and derive 
in the limit the following stochastic Hele-Shaw problem 

' Av = 0 in V\T(t), t > 0, 
d n v = 0 on dT>, 
v = XH + W on T(i), 


(2.5) 


V = \{d n v + -d n v ) on F(t), 


r(0) = r 0 , 

where again H and V are the mean curvature and velocity respectively of the zero level surface 
T(t). For positive e > 0 the domain T> admits the following disjoint decomposition 

v = vt(t)uv~(t)uv I E (t), 

where 

u e ps 1 for xET>f(t) and u e « — 1 for x G V~ (t). 

Moreover, T>l{t) is a narrow interfacial region around T(t) with thickness of order e where u £ is 
neither close to +1 nor —1. 

In particular, we construct an inner solution close to the interface, and an outer solution away 
from it. Using the appropriate matching in orders of e, we formally pass to the limit and derive 
the corresponding free boundary problem. To avoid additional technicalities, we also assume that 
the interface T does not intersect the boundary. In terms of simplicity of notation, we drop the 
subscript e in all the calculations that follow. 

2.2. Outer expansion. We consider that the inner interface is known, and seek the outer expan¬ 
sion far from it, i.e., an expansion in the form 

u = Uq + eu\ -\ -, 

v = Vq + £V\ H-, 

where ‘+...’ denote higher order terms and uq,u\, ... ,vq,vi, ... are smooth functions. We insert 


the outer expansion into the second equation of the stochastic system 


and obtain 


1 


Uo + £V\ + 0{e ) — — —{f'(uo) + Ef"('Uo)(ui + W) + (D(e )) 

+ £A('Uo + £U\ + 0[e 2 )) + e 2 AW + 0(£ 2 ). 
1 


( 2 . 6 ) 


First collecting the terms of order J in (12.61) . we arrive at 

f\u 0 ) = 0. 

Thus, we get as in Remark 4.1, (1) of (2] 

uo = ±1 • 

In the second step, we collect the 0(l)-terms in (|2.6D and derive 

u 0 = -f"(uo){ui + W) + 0(£ 2 ). 

We plug now the outer expansion into the first equation of (12.41) , and obtain 


dt(uo + £Ul + C>(£ 2 )) = -A(u 0 + £V\ + 0(£ 2 )) . 
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As uq is a constant in the outer expansion we have dtuo = 0, and thus, collecting the 0(1) terms 
yields 

—Avo = 0 . 

Collecting finally in the third step all (9(e)-terms we arrive at 

d t u\ = -Am. 


2.3. Inner expansion. Let x be a point in V that at time t is near the interface T(f). Let us 
introduce the stretched normal distance to the interface, z := where d(x,t) is the signed distance 
from the point x in T> to the interface T(i), such that d(x,t ) > 0 in T>f and d(x,t ) < 0 in T>~. 
Obviously T has the representation 


T(f) = {x € V : d(x, t ) = 0}. 


If r is smooth, then d is well defined and smooth near T, and |Vd| = 1 in a neighborhood of T. 
Following [2] and m, we seek for an inner expansion valid for x near T of the form 


u = q 


x , 7 ) +e q(^,x, <) + .... 


d(x, t) 


v 


q 


d(x, t ) 
£ 



+ £Q 



+ • • • , 


where again ‘H-’ denote higher order terms that we neglect and q,Q,... ,q,Q,... are sufficiently 

smooth. It will be convenient to require that the quantities depending on z ,x, t are defined for x 
in a full neighborhood of T but do not change when x varies normal to T with z held fixed, [47]. 
We insert the inner expansion into (12.4jL . utilize that |Vd| 2 = 1, and thus obtain the following 
expression 

(2.7) q + £Q + 0(£ 2 ) = --(f'(q) + £f"(q)(Q + W))+£(^Ad+^- + d z QAd+^-)+£ 2 AW. 

We collect the terms of order an d derive 

dzzq ~ f'(q) = 0 . 


By matching now the terms of order 0(1) in (12.71) . we obtain 


q = —f"(q)Q + d zz Q - f"(q)W + d z qAd, 


or equivalently 

(2.8) q - d z qAd = d zz Q - f"(q)Q - f"(q)W ■ 

We define the linearized Allen-Cahn operator 

CQ = d zz Q - f"(q)Q . 


Then (12.81) is written as 

(2.9) q-d z qAd = CQ- f"(q)W . 

This equation is solvable if for any % £ Ker(£*) it holds that x -L (q 
equivalently if 


(2.10) 


X -(q-d z qAd + f"(q)W)dz = 0. 


d z qAd + f"(q)W ), or 
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Obviously, for any x on T it holds that d(x,t ) = 0 and A d(x,t) = H(x,t). Replacing in (|2.10p we 
obtain the following sufficient condition on the interface F: 

(2.11) q = XH-f"{q)W. 

Plugging the inner expansion into (I2.4p i we obtain 

( 2 . 12 ) ^l dt + dz Qd t + 0{e) = -(^Ad+^- + d z QAd+^-). 

We collect the terms of order (D(e~ 2 ) and arrive at 

<9 zzQ = 0, 


which implies that for some functions a and b 

q = a(x, t)z + b{x, t ). 

To proceed further, the matching condition for the inner and outer expansions must be developed. 
In general, these are obtained by the following procedure (see m)- Fixing i 6 T, we seek to match 
the expansions by requiring formally for z —> oo 

q + £Q + 0(e 2 ) = uo + + 0(e 2 ), 

and thus in order 0(1) 

vo = lim q = lim ( a(x,t)z + b(x,t)). 

z—>oo z—>-oo 

We obtain a = 0 and thus, q = b. Hence, utilizing (|2.11j) we have that on the interface 

v 0 = XH- f"(q)W = XH + W , 

where we used that q solves the Euler-Lagrange equation 

- + f'(q(z)) =0, Rl, 

lim q{z) = ±1, ?(0) = 0, 

£—>■±00 


while on the inner interface with 2 = d/e = 0 we have f"(q) = 3 q 2 — 1 = —1 since g(0) = 0. 

What is still missing is the evolution law, which should come from the inner expansion. From 
(12.121) . we collect the terms of order 0(l/e) and obtain 

d z qd t = -d z qAd - d ZZ Q . 


Recall that —dt = V, while Ad = H (see for example 0) and integrate over z from —00 to 00 to 
derive 

/*oo /*oo 

- / d z qV dz = — d zz Q dz . 


From the matching conditions we get 

g(+oo) = 1 and q{— 00 ) = —1 . 


Hence, we have 

V = \[d z Q{+o o) - d z Q{- 00 )] . 

Thus the stochastic Hele-Shaw problem (12.51) is established formally as the sharp interface limit, in 
the case <7 = 1. 
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Remark 2.1. In the case a > 1, we follow the same construction of inner and outer solutions as 
above and obtain in the limit, the deterministic Hele-Shaw problem, m) 

' Av = 0 in V\T(t), t > 0, 
d n v = 0 on dV, 

< v = A H on T(t), 

V = d n v + - d n v~) on T(t), 

. r(o) = r 0 , 

where H and V are the mean curvature and velocity respectively of the zero level surface r(i) 
contained in the interfacial region T>l(t). 


Remark 2.2. Note that the change of variables 

u e := u e - s a W, 


implies that 


u £ := u e + 0(e a ). 


Thus, only if a > 1, u is permitted to be expanded as in the presented inner expansion using q. The 
key difference to the deterministic analysis is that in the nonlinearity we have jf(u e + e ff W). In 
case a > 1 there is no contribution of IN in an asymptotic expansion in terms of order 0(1) and 
0(1 /e), while for a = 1 there is an impact of IN on terms of order 0(1). 


Remark 2.3. When 0 < cr < 1 the strategy presented in this section fails. For this case, we might 
think of avoiding the change of variables and apply the formal asymptotics of [TO] directly instead. 
If we split the noise term W in m, we obtain the following stochastic system with two noise 
sources (recall that e CT VV = AV" + W) 

' d t u = -A v + W, 

v = -^M + £ a u + V, 

£ 

The sharp interface limit should coincide to (12.31) , but for lim G\, lim G 2 replaced by 

£->■0+ e->0+ 

lim and limy(-,e), 

£->0+ £^-0+ 

respectively. When 0 < a < 1 we would obtain that both these limits are 0 and thus the limiting 
problem is a deterministic Hele-Shaw problem without the contribution of the noise. But this is a 
very dangerous reasoning, as the noise terms V and W, even if they would be e-independent are 
not of order 0(1), and we would still expect an impact of the noise terms on the limiting problem. 


3. The sharp interface limit 

The main result of this paper is that dm as e tends to zero, may have a deterministic or a 
stochastic profile depending on the strength of the additive noise in terms of e. Only large noise 
perturbations with a = 1 generate a stochastic limit problem. Here, we discuss the limit for smaller 
noise strength. 
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Let us first precisely state our problem. We assume that noise is induced by the formal derivative 
of a Q-Wiener process W in a Fourier series representation (see my, for simplicity, the only e- 
dependence will appear in the noise strength, and thus, for the rest of this paper we shall use the 
notation dW{x,t) for the additive noise, where a £ R is a scaling parameter. 


Assumption 3.1. Let W be a Q-Wiener process on some probability space (0, A, P) such that 

W(t) = £ a k /3 k (t)e k , 
fce N 

for an orthonormal basis {e k ) k ^ in L 2 (V), independent real-valued Brownian motions (/3 k )keN on 
(fi, A, P), and real-valued coefficients a k such that Qe k = a k e k . Furthermore, we assume that the 
noise has some weak smoothness in space, i.e., Q satisfies 

(3.1) trace(A _1 Q) < oo. 

To deal with a mass-conserving stochastic problem, we impose the condition 



W(t)dx = 0 . 


Note that (13.111 implies that the Wiener-process W(f) is lL _1 ((D)-valued. This is the minimal 
requirement for the approximation theorem presented in the sequel; we might need more regularity, 
in order to have the stochastic Hele-Shaw limit problem well defined, or while performing the formal 
asymptotics. 

Recall m in Ito-formulation 

(3.2) du £ = A(—eA u £ + e _1 f' (u £ ))dt + £ c dW{x , t), 


associated to Neumann conditions on the boundary of V, so that the equation is still mass conser¬ 
vative. 

The following theorem for the existence of solutions is well known. See for example m- 

Theorem 3.2. Let T> be a rectangle in dimensions 1,2,3. If Q = I or trace(A~ 1+s Q) < oo for 
5 > 0, then the following holds true: 

(1) if uq is in H^ifD), there exists a unique solution for the problem (11.111 in C([0,T]; iL~ 1 ((D)), 

(2) if uq is in L 2 (V), then the solution for the problem (11.111 is in L°°(0, T; L 2 (V)). 


Note that the previous theorem could be extended for general Lipschitz domains in dimensions 
2 and 3 under some additional assumptions of minimum eigenfunctions growth, cf. the arguments 
in [9]. For the analysis underlying our results and to avoid technicalities, we will for the remainder 
of the paper always assume: 

Assumption 3.3. For any initial condition in H~ l (fD) there exists a unique solution for the 
problem mm in (7([0,T]; H l (fD)), which is sufficiently regular such that we can apply ltd-formula 
to the H^ 1 -norm. 


The regularity required for this assumption is usually straightforward to verify. For example m 
applies the Ito-formula to a finite dimensional spectral Galerkin-approximation and then passes 
to the limit. Introducing the chemical potential v £ , the equation is as in the formal derivation 
rewritten as a stochastic system. Indeed, 

for T > 0 let Vt := V x (0, T), 
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then (13.21) is written as 


(3.3) 


du £ = — A v £ dt + £ a d\V 
v £ = - \ f'{u £ ) + eAu £ 


subject to Neumann boundary conditions 


in T>t, 
in T>t, 


du £ 

dn 


dA u £ 
dn 


= 0 


on dV x (0, T ) 


Our main analytic theorem considers a sufficiently small noise resulting to a deterministic sharp 
interface limiting behavior. In particular, we analyze the case 


(T » fj 0 = 1, 

where ctq is the conjectured borderline case, where according to our formal calculation the noise 
has an impact on the limiting model. Under some assumptions on the initial condition w e (0), the 
limit of u £ and v £ as e —> 0 solves the deterministic Hele-Shaw problem on a given time interval 
[0,T], We will state the precise formulation of this argument in Theorem 13.101 and then present 
the rigorous proof. 

While our rigorous result can only treat very small noise, the formal derivation of Section 12.11 
motivates the following conjecture implying a stochastic sharp interface limit: 


Conjecture 3.4. For a = 1 the limit of u £ and v £ solves the stochastic Hele-Shaw problem 

Av = 0 in V\T(t), t > 0, 
d n v = 0 on dV , 
v = XH + W on r(t), 

V = d n v + - d n v ~) on T(f), 

r(o) = r 0 . 

Remark 3.5. In Section 2 by formal asymptotics, we only presented an indication for the correct¬ 
ness of the conjecture. A rigorous proof of this conjecture remains open at the moment. We hope 
to attack the problem to its full generality in the near future. 

Remark 3.6. Note that W is a Wiener process, and the equation v = A H + W on T(t) appearing in 
has a rigorous mathematical meaning in terms of functions. In fact, no noise is present, while 
a random equation appears on V\T(t) in the following sense. For any given t, T(f) is defined by its 
velocity V and thus is known, and the unknown function v on T(t) is a stochastic process. Thus, 
the problem for fixed t is posed in between the inner boundary T = T (f) and the outer boundary dV 
as follows 

(3.5) Av = 0 in T>\T, d n v = 0 on dV, v = A H + W on T, 

the inner boundary condition u|r being a realization of a t-dependent stochastic process. 

3.1. Statement of the Main Theorem. In this section, we shall state the main analytic theorem 
of this paper, concerning the sharp interface limiting profile for sufficiently small noise strength. To 
approximate the stochastic solution we use the same approximations uf and vf as in [2] proposed 
in the absence of noise. For a precise definition see further below. In our proof we follow the ideas 
of the proof of their Theorem 2.1, and need to adapt the analysis to the presence of noise. 


(3.4) 
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The main difference concerns the noise in the equation for the residual 
(3.6) R := u £ — uf. 

We show bounds for this error in our main Theorem 13.101 below, while in Remark l3.11l we comment 
on the smallness of the noise necessary for the main result. Unfortunately, we are quite far away 
from treating the case a> 1 close to 1. Another problem is the weak estimate on the spectral 
stability of the linearized operator stated in Proposition 13.91 See also Remark 13.131 


Assumption 3.7. Let the family {r(t)}ig[o,T] °f smooth closed hypersurfaces together with the 
functions {r(t)}te[o,T] be a solution of the deterministic Hele-Shaw problem (i.e., equation (13.41) 
with W = such that the interfaces do not intersect with the boundary dT>, i.e., T(f) C V for all 
te[0,T}. 


With T from Assumption 13.71 the authors in [2] construct a pair of approximate solutions 
(uf, vf ), so that T(t) is the zero level set of uf(t ), which satisfies 


(3.7) 

for boundary conditions 


duf = — Avfdt 


in T>t j 


= - 7 f\ u e ) + eAu e + r, 


in T>t, 


dut 


dA uf 
£ =0 


on ()'D. 


dn dn 

We recall that uf approximates the deterministic version of equation (11.11) (i.e., for W = 0). The 
error term rf is bounded in terms of e, and depending on the smoothness of T and the number of 
approximation steps, the bound on rf can be arbitrarily small. For details see relation (4.30) and 
Theorem 4.12 in [2]. 

We will summarize the results of [2] that we need for our proof in the following Theorem. 


Theorem 3.8. Under the Assumption \3 . 7\ for any K > 0 there exists a pair (uf, vf) of solutions 
to ( 13 . 71 ) . such that 

\\ r t\\c°{v T ) < Ce k 2 . 

Moreover, it holds that 

W'v'e ~ v \\c°(p T ) < Ce , 

and finally for x away from T(t) (i.e., d(x,T(t)) >ce) 

\uf(t, x) — 1| < Ce or \uf(t,x) + \\<Ce. 


We present now the following spectral estimate, useful in our proof; we refer to f22] for dimensions 
larger than two, and to d] for dimension two. Unfortunately, this estimate is also the key problem 
to extend the approximation result beyond time-scales of order 1. 

Proposition 3.9 (Proposition 3.1 of [2j). Letuf be the approximation given in Theorem \3.8\. Then 
for all w G i7 1 ((D) satisfying Neumann boundary conditions such that f v wdx = 0, the following 
estimate is valid 

I [e|Vu!| 2 + ~ f (uf)w 2 }dx > — Co||Vu;|| 2 2. 
j v £ 

Our main theorem provides bounds for the residual R: 
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Theorem 3.10. (Main Theorem) Let Assumption \3. 1\ for the noise a,nd \3.3\ for the existence of 
solutions be true. 

Fix a time T > 0 and any pG (2,3], a radius e 7 with 

' 2p + d(p-2) p + 2' 

2 p — d(p — 2) p . 

and a noise strength e a with 

2 p + dip — 2 ) p + 2 
a> ^ + 2p-d(p-2) - — - 

Then for some small k > 0 and a generic constant C > 0 there is for all large £ > 0 a constant 
Ci such that the following is true: 

For all families {F(t)} te [ 0 ,'Z’l of closed hypersurfaces satisfying Assumption \3. r !\ with correspond¬ 
ing approximation uf and vf defined above in Theorem I3.<gl and u £ solution of the stochastic 
Cahn-Hilliard equation (13.21) such that it £ (0) = uf( 0), the following probability estimates hold: 

^ (II-R||lp([0,T]x 23) < e 7 ) > 1 ~ Ci£ e , 

IP (ll-RHico^r.fr-i) < + e"-*-*]) > 1 - C t e* , 

and 

If (ll«lli>(0,T,ff.) < C[c-'- 2 'r+*r + J-1+-+7-]) > 1 _ c,e ‘, 

where R := u e — uf is the error defined in m- 

Remark 3.11. Let us remark that in dimension d = 2 one can easily check that we obtain the 
smallest possible value both for a and 7 for p = 3. In that case 7 > 6 and a > 23/3. This is in 
well agreement with the 7 derived in [ 2 ], but unfortunately we can only consider very small noise 
strength. But it seems that using the H~ 1 -norm and spectral information available there is no 
improvement possible. 

For dimension d = 3 again the noise strength is small, but the result is not that clear. While the 
smallest value for 7 is still attained at p = 3 (with 7 > 6 and a > 11) we obtain the smallest value 
of a for some p < 3. 

Remark 3.12. Let us remark on the fact that we take R( 0) = 0. We could take more general 
initial conditions u £ { 0) for the Cahn-Hilliard equation. Looking closely into the proof of Theorem 
\3.10\ we could allow for an initial error R( 0) with ||R(0 )||^_ 1 = 0(e P7_1 +e°‘ +7 ~ K ). 

Remark 3.13. Let us state two main problems with the approach presented. 

First, the spectral estimate in Theorem \3.9\ yields an unstable eigenvalue of order 0(1). This 
immediately restricts any approximation result to time scales of order 0(1). But we strongly believe 
that this eigenvalue represents only a motion of the interfaces itself. One would need spectral 
information orthogonal to the space of all possible approximations uf, which are parametrized by 
the hypersurfaces T. But this does not seem to be available at the moment. 

Moreover, later in the closure of the estimate we can only allow a > oq large enough, i.e., for 
sufficiently small noise strength. Here, an additional problem is that the H^ 1 -norm is not strong 
enough to control the nonlinearity, and from the spectral theorem \3.9[ we do not get any higher 
order norms that would help in the estimate. Nevertheless, if we start with higher order norms like 
L 2 , for instance, then there are no spectral estimates available at all. 
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4. The proof of the Main Theorem 13.101 


4.1. Idea of Proof. For the proof we define for p G (2,3] and cr > 7 > 0 (both fixed later), the 
stopping time 

(4.1) T e := inf |f G [0, T] : ( j\\R{s)\\ p LP ds) 1/P > 

where the convention is that T e = T if the condition is never true. 

The general strategy for the proof of the main theorem is the following: 

(1) Use Ito-formula for d\\R\\ 2 H _ 1 . 

(2) Consider all estimates up to T e only. 

(3) Bound the stochastic integrals (at least on a set with high probability). 

(4) Show that T e = T with high probability using the bound derived for f* up to T e . 


4.2. A differential equation for the error. Let us first derive an SPDE for R from (13.6(1 . using 
(13.71) and (13.3() . as follows 

dR = du e — duf = A vfdt — A v e dt + e a dW 


(4.2) 


- - e M\uf) + eA 2 u? + Arf + i A/'(u e ) - eA 2 u, 


1 


Af\uf + R)-Af\u-£) 


dt T 


- eA 2 R + Art 


dt + £ a dW 
dt + E a dW. 


4.3. The H '-norm of R. The approximate solutions uf and vf are, by their construction, 
functions in C 2 (T>t ), while uf satisfies for all t G [ 0 ,T] 


/ uf(t)dx = 0 . 

Jv 

Since (13.21) is mass conservative, we can conclude that mass conservation also holds for R, i.e., for 
all t G [0, T] 

/ R(t)dx = 0 . 

Jv 

Observe that the operator —A is a symmetric positive operator on the space 


R 2 := ju; G C 2 {V) : j w dx = 0 and d n w = 0 on dnj. 

Therefore, by elliptic regularity, the operator —A : R 2 —»• L 2 is bijective. So, we can invert it and 
for any t G [0, T] there exists a unique ij)(t) G R 2 such that 

(4.3) — A x[)(t) = R(t), or equivalently (— A)~ 1 R(t.) = ij(t). 

With the scalar product (■, ■) in L 2 the U -1 -norm of R is given by 

\\R\\h-i = ||(-A)- 1 / 2 i?||i 2 = ||(-A) 1 /^||2 2 = \\ vn i 2 = { ^ R) . 


Since 

{dJj,R) = (-A dR,R) = (dR,~AR), 
considering the Ito-differential, we obtain 

y\\R\\ 2 H -i = {'ip, dR) + i {dip, dR) = {ip, dR) + ie 2a ((-A )~ l dW, dW) 
= {iP,dR) + ie 2,7 tr(Q 1 /2(_ A) -iQi/2 ) _ 


(4.4) 
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Here, by Assumption 13.11 the trace in the previous estimate is bounded. So, using (13.31) and (13.7p . 
we arrive at 


(4.5) 


(ip, dR) = (ip, d(u £ - uf)) = (ip, (~A)(v £ - vf)dt + e a dW) 
= {R , ( v £ — vf))dt + £ a (ip, dW). 


Using again (13.311 and (13.711 in order to replace the v’s, yields the following equality 

(ip, dR) = — e~ 1 (R, f'(u e ) — f'(uf))dt + e(R, A (u £ — uf))dt — (R, rf)dt + £ a (ip, dW). 
= -e _1 (R, (f(u £ ) - f(uf)))dt - £\\VR\pdt - (R,rf)dt + £ a (ip,d\V). 


(4.6) 


Definition 4.1. For any positive integer p, we define the L p -norms 



Also, we denote by || • || the usual L 2 (ID)-norm and by || • \\u> the L V (T>)- norm. 

Applying Taylor’s formula to expand f(u £ ) around uf, with residual J\f(uf, R), we have 

f(u £ ) - f(uf) = f"(uj)R+N(utR) ■ 

The crucial bound for the nonlinearity in the residual is the following result from Lemma 2.2 of [2|. 
It is based on a direct representation of the remainder M in the Taylor expansion together with 
the fact that uf is uniformly bounded. 

Lemma 4.2. Let p £ (2,3] and q such that ^ | = 1, then it holds that 

(4.7) -J £- l U(utR)R<c£- 1 \\R\\l v . 

Thus, we obtain 

-- <K. (/'(«<) -/'W))> = -\ rj "( u a )r) _ i 

(4.8) e £ £ 

<- ; {iJ,/"(^)/J)+ce- 1 ml,v 

Relations (14.411 . (14.6p and (14.8H yield the following first key estimate 

(4 g \d\\S7iP\\ 2 + e\\VR\\ 2 dt + - e (R, f"(uf)R)dt 

< c£ _1 ||R||p iX ,df + \\R\\ p ,v\\r £ \\ q ,vdt + £ a (ip, dW) +C tr £ 2a dt. 

From this a-priori estimate, we now can derive a uniform bound for || V'0|| and later a mean square 
bound on ||Vi?||. Both estimates still involve the L p -norm of R on the right hand side, and we use 
the stopping time T £ to control this. 


4.4. Technical Lemmas. We first need the following Lemma of Burkholder-Davis-Gundy type 
for stochastic integrals. Recall the stopping time T e from ( 14 . 111 . 

Lemma 4.3. Let f be a continuous real valued function, and A ip = R as before. Then for all 
k > 0, £ > 1 there exists a constant C = C(£, T, k ) such that 

pf sup | f f(ip,dW) | > e 7_K 
v [o , T e ] Jo 


< Ce^\\ft 2pnp _ 2) 
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Proof. We shall use the Chebychev’s inequality. Thus, we need to bound the moments first. Apply¬ 
ing Burkholder-Davis-Gundy inequality (using that A _1 QA _1 is a bounded operator by assump¬ 
tion), we obtain 


E sup 

[0,T e ] 


f(s){'ip(s),d'W(s)) 


<CgE 

<CE 

<CE 


I f 2 {s)('ip(s),Q'ip(s))ds 
o 

rTe a /2 


1/2 


>0 


f (s)\\R(s)\\hds 


\\R(s)\\ p LP ds 


2/p 


1/2 


f2p/(p-2) 


(p-2)t/(2p) 


<C||/||' 1i/( „_ 1)£ A 

Here, we applied Holder’s inequality and the definition of T e < T. 

Furthermore, using Chebychev’s inequality, we obtain the result as follows 

P( sup | f f(if,dW )| > e^) < sup f f(if,dW) * < Ce £K \\f\f L2p/{p _ 2) . 

V [ 0 ,Te] Jo 7 [0 ,T e ] Jo 


□ 


Now we present the following stochastic version of Gronwall’s Lemma. 

Lemma 4.4. Let X, Pi, A be real valued processes, and Q be a Hilbert-space valued one. Further¬ 
more, assume that 

dX := [XX + P 1 \dt + (G, dW), 

and that 

F\ < -F 2 ■ 

Then the following inequality holds true 

X(t)<e m X( 0)+ / e A(t) " A(s) 7' 2 ( s )d. s+ f e Ait) ~ A{s) (G(-s), dW(s)), 

Jo Jo 


for 


Aft) := f A (s)ds. 
Jo 


Proof. We define 

Y(t) := X(t)e~ A{t l 

By the definition of the process Y, we obtain easily 

dY = e~ A dX - A Ydt = e" A J'i dt + e~ A (G , dW), 

and 

Y{t) = Y{ 0)+ fe-^Fi [s)ds+ f e~ A ^{G(s), dW{s)) 

Jo Jo 

<A(0)+ [ e~ A ^F 2 (s)ds+ [ e~ A ^{G(s),dW(s)). 

Jo Jo 

Multiplying the inequality with e A J\ and using the definition of Y, we derive the stated stochastic 
version of Gronwall’s inequality. □ 
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4.5. Uniform Bound on 'S/ip. Using the spectral estimate of Proposition 13.91 we get from (|4.9I) 
(4.10) d\\ViP\\ 2 < [c||V^|| 2 + ce- 1 ||^||^ + 2||ii|| PiC ||r £ A ||^ + Ce 2CT ]dt + 2 e ff (i/;,dW). 

We apply Lemma fOl Since i?(0) = 0 implies Vip(0) = 0, this yields 


l | Vi /’(*)|| 2 < f 

Jo 


r,C(t—S) 


08 ll-^llp,X) + \\ R \\ P ,A r e \\ q,V + Ct T e 


< e 


•+/ 

f C£~ 1 \\R\\ P V + \\R\\ p ,v\\rt\\ q ,v + C tY £ 2a ds + e a e CT [ 
Jo 1 J Jo 


s C(t—s) a 


£ a (ip, dW(s)) 


e Cs (ip,dW(s)) 


Furthermore, from Lemma 14.31 we obtain on a subset with high probability 

pt 


sup 

tE[0,T e 


e Cs (ip(s),dW(s)) 
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< Ce 1 ' 


Thus, we arrive at 

l|V'0(t)|| 2 < Ce~ l \\R\\ p pVt + C\\R\\ p , Vt \\r%v t +Ce 2 °t + Ce° + ^ 
+ £ \\r e || q , V T + £ ~r £ 


(4.11) <C[£ P1 ~ l -^^\\r- A \\ ~ -Lr 2ff J-c ff +7-«l 


< C^- 1 + £ a+ry - K ], 

where we used that 7 < <7 and that k is sufficiently small, together with Theorem 
verified the following Lemma: 


Finally, we 


Lemma 4.5. For all p € [2,3), a > 1 , k > 0, and 7 < cr we have 
(4.12) <C [£^- 1 + ^-«], 

with probability larger that 1 — C^ for all l > 0. 

4.6. Mean Square Bound on VR. We return to relation (14.91) and shall use the estimate 

.-1 


—e 


0 Jv 


f\u*)R 2 dx<£~Vv\\R\\l Vv 


presented in j2] on p. 171. Its proof is based on Holder inequality together with the fact that the 
set where the value of uf is not close to either +1 or —1, has a small measure. More precisely, the 
measure is controlled by 


measured 


j{(M) e < oj < Cer, e € (0,1]. 

Therefore, integrating ()4.9|) and using Vip( 0) = 0 (since R( 0) = 0), we arrive at 

4VR\\iv t < + \\R\\ P ,v t \\r%,v t + £ a f(ip,dW) + C^t. 

Jo 

Revoking again Lemma 14.31 we obtain on a set of high probability 

ellVRH!,^ < e - 2/ '||fl|g,^ + c£- 1 ||R||^ + || R\\ p ,v Te \V%,v T + + C tr £ 2 °T, 

where we used that T e < T. Moreover, the definition of T e implies for all t £ [0, T e \ 
ellVRUl^ < £~ 2 ' p £ 2 ^ + C£~ l £^ + ^WrfW^Vr + £ a+1 ~ K + C£ 2 F 
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_er+7—K1 


Here, the constant depends on the final time T. Using again 7 < a and k sufficiently small, together 
with Theorem 13.81 we obtain 

(4.13) e||Vi2||l )2>re < C[e~ 2 ^e 2 ^ + £"V 7 

Note that as p > 2, a short calculation shows that 

r — 2 /p 27 


£ > £ _1 £ P7 


1 

- < 7, 
P 


which we assume from now on, as we expect both 7 and a to be bigger that 1. We verified the 
following Lemma: 

Lemma 4.6. For all p £ [2,3), k > 0, a > 1, and 7 £ (|, a) we have 

(4.14) l|fl(t)lli. ( o,T„H>) < C[e -I-2/,p+27 + e - l + a +r- K j 

with probability larger that 1 — Cfor all £ > 0 . 


4.7. Final step. In the final part of the proof it remains to show that T e = T on our set of high 
probability. Thus, we shall use our estimates of the previous two Lemmas to show that ||i2|| Pi x> is 
not larger that £ 7 . 

Observe first, that the following trivial interpolation inequality holds true 

(4.15) \\R\\ 2 2 v = ~ [ RAifdx= [ 'VKVifdx < ||V-R|| 2 ,x>||V'0||2,x>. 

Jv Jv 

We use the Sobolev’s embedding of H a into L p with a := d(^ — 4) = 2 p , an d then interpolate 
H a between L 2 and H 1 . We need a £ [0,1], which is assured by 2 < p < 3 < ^ 2 ) • This gives, 

\\R\\ p ,v<C\\R\\ h ° <C'||i?||^||Vi?||^. 

Thus, using (14.151) we obtain 

2p—d(p—2) d(p— 2 ) 

\\R\\ P P ,v<C\\R\\^^\\VR\\^r 

2p—d(p—2) d(p— 2 ) . 2p—d(p— 2 ) 2p—d(p—2) d(p-2)-\-2p 

< C'||V'0II^^ l|vi?||^ 3 = c'livV’ll^T 2 . 

Integration yields 

2 p— d(jp— 2 ) d(p— 2 )+ 2 p 

PIIU < Cl sup IIVV’II^T* • IIV^II^* . 

[04] 

Now we use (14.121) and (14.141) and arrive at 


(4.16) 


\\R\\lv < ci 


£ P7-1 _|_ £ cr+7-K 


2 p-d(p- 2 ) 


1 — 2/p+27 _|_ £ — l+fT+7 — K 


d(p— 2 )+ 2 p 
8 


Thus, pulling out £ 27 from both brackets, and as 7 — <7 > 0 and n is small, we arrive at 


e^\\R\\p,v<C- 

<c- 


r(P~ 2)7-1 _)_ e <T-7-K 


2 p—d(p— 2 ) 


-l-2/p _|_ £ -l+ 0 — 7 -K 


d(p— 2 )+ 2 p 


r(P—2)7— 1 _)_ £ <T-7-K 


2 p— d(p— 2 ) 


^_ £ _(l +2 /p )M±2p 
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By Lemmas 14.51 and 14.61 the previous bound holds with probability larger than 1 — Cisr. In order 
to show that T e = T holds with high probability, we need to prove that the right hand side of the 
previous equation is smaller than 1. 

As the second factor is larger than one, we need the first one to be smaller than 1 and to 
compensate the larger factor. Hence, we need 


7 > 


p — 2 . 


1 + 


2 p + d{p - 2) 
2 p — d{p — 2) 


p + 2i 1 

- - > - 

p J p 


and provided k is sufficiently small 


a > 7 + 


2 p + d{p — 2) 
2 p — d{p — 2) 


p + 2 

->7- 

p 


Hence, the proof of theorem is now complete. 
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